Magnetic Field Effect on the Pseudogap Temperature within Precursor 

Superconductivity 
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We determine the magnetic field dependence of the pseudogap closing temperature T* within a 
precursor superconductivity scenario. Detailed calculations with an anisotropic attractive Hubbard 
model account for a recently determined experimental relation in BSCCO between the pseudogap 
closing field and the pseudogap temperature at zero field, as well as for the weak initial dependence 
of T* at low fields. Our results indicate that the available experimental data are fully compatible 
with a superconducting origin of the pseudogap in cuprate superconductors. 



PACS numbers: 74.20.-z, 74.25-q, 74. 25. Ha 

High-temperature superconductors are characterized 
by a pseudogap phase, whereby a temperature T* for 
the closing of the pseudogap is identified in addition to 
the (lower) superconducting critical temperature Tc- The 
origin of this pseudogap phase remains controversial and 
has mainly been ascribed to two alternative scenarios, 
namely, to the occurrence of a "competing order param- 
eter" (due, e.g., to an underlying antiferromagnetic or 
charge order) or to the presence of "precursor supercon- 
ductivity" (due to strong pairing attraction above Tc). 
The issue is thus to decide whether pseudogap and super- 
conductivity (better, precursor pairing) represent related 
or different phenomena. 

In this context, magnetic effects are good candidates 
for distinguishing between the two scenarios. It turns out 
that T* and Tc behave differently as functions of hole 
concentration and magnetic field H, in the sense that, 
in the underdoped region, T* is weakly dependent on H 
while Tc shows more pronounced dependenceo. A simi- 
lar behavior is observed at optimal dopingnl, while in the 
overdoped_region the pseudogap behavior is strongly field 
dependentcl. The above distinction between the two tem- 
perature (or energy) scales has sometimes been quoted 
as evidence that the two scales have different physical 
originQiJ. 

Recently, by applying magnetic fields up to 60 T, a 
systematic determination of the pseudogap closing field 
Hpg (above which the pseudogap phase is destroyed at 
any temperature) has been made using data for several 
doping valuesEi, showing strikingly different dependen- 
cies from the characteristic fields of the suaerconduct- 
ing state. In addition, it has been foundO that Hpg 
and T* at zero field are related (within experimental 
uncertainty) through the simple Zeeman-like expression 
g^sHpg = ksT* , where g = 2 is the electronic g-factor, 
(jiB the Bohr magneton, and ks the Boltzmann constant. 
This expression suggests that the magnetic field couples 
to the pseudogap by the Zeeman energy of the spin de- 
grees of freedom, thus entailing a predominant role of the 



spins over the orbital frustration effects in the formation 
of the pseudogap. 

Within the framework of the pairing scenario, a cal- 
culation of the effect of a magnetic field on T* and Tc 
has already been presented by including orbital mag- 
netic effects only but no Zeeman coupling in a continuum 
modeD, with the conclusion that orbital effects should 
contribute in an important way to the field dependence 
of T*. 

In this Letter, we calculate the magnetic field depen- 
dence of T* within a precursor superconductivity sce- 
nario, by taking explicit account of the Zeeman split- 
ting for the spin degrees of freedom, and provide spe- 
cific comparison with the data of Ref. |^. To this end, 
we initially consider the Zeeman splitting only, within 
an attractive Hubbard model for an anisotropic (lay- 
ered) structure with a d-wave pseudogap. At a sec- 
ond stage, we include also orbital effects in a continuum 
model, to assess the relative importance of Zeeman and 
orbital couplings. Our results are that, when pseudo- 
gap effects are weak (i.e., at low pairing interaction), T* 
depends markedly on H, while at stronger pairing T* 
is weakly field dependent. In addition, at the doping 
levels relevant to experiments, we reproduce quantita- 
tively the expression gfisHpg = ksT* proposed by the 
authors of Ref. ||. Specifically, we find that: (i) Inclusion 
of the sole Zeeman splitting for an attractive Hubbard 
model with an anisotropic (layered) structure recovers 
the above expression within about a factor of two (which 
we attribute to the difference between the temperature 
T* as identified within precursor superconductivity and 
the energy scale associated with pair breaking); (ii) The 
contribution of the orbital effect (estimated for a con- 
tinuum model) does not significantly affect the results 
obtained for Hpg with the sole inclusion of the Zeeman 
effect, in the intermediate-coupling regime relevant to ex- 
periments; (iii) In the small H region, T*{H) is weakly 
dependent on H and somewhat sensitive to the inclusion 
of the orbital effect. Our results imply that the magnetic 
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field effects in the pseudogap phase, as determined by the 
experiments reported in Ref. ||as weU as in Refs. and 
^, are fuUy compatible with precursor superconductivity 
as the origin of T* . 

The relative importance of the Zeeman and orbital ef- 
fects on ifpg could be assessed by applying a magnetic 
field, alternatively, parallel and orthogonal to the lay- 
ered structure of the high-temperature superconducting 
materials, thus suppressing the orbital contribution when 
the magnetic field is parallel to the layers. Although no 
experimental indication in this sense yet exists, the rela- 
tion between iJpg and T* found in Ref. || constitutes per 
se an empirical suggestion that the Zeeman effect domi- 
nates over the orbital frustration, at least as far as iJpg 
is concerned. 

For the above reasons, we consider initially the Zeeman 
effect only, which can be evaluated for arbitrary magnetic 
field and for which lattice effects (even with anysotropy) 
and the d-wave character of precursor pairing can be 
readily included. To this end, we consider an attractive 
extended Hubbard model, with nearest-neighbor (t||) and 
next-nearest neighbor in-plane hoppings, nearest- 
neighbor (t±) out-of-plane hopping, and an in-plane at- 
traction (with strength V < 0) between nearest-neighbor 
sites. This model is believed to capture the essential 
physics of cuprate superconductors in the nearly opti- 
mally doped region. 

In the absence of magnetic field, a pair breaking tem- 
perature is obtained by supplementing the standard BCS 
gap equation with the particle-density equation to deter- 
mine the chemical potentialEI, since its value is no longer 
pinned at the Fermi level as soon as the attraction is suf- 
ficiently strong. When this happens, the above temper- 
ature is no longer related to the superconducting critical 
temperature, but rather signals the formation of (quasi) 
bound pairs. In this context, the pair breaking temper- 
ature is then associated with the crossover temperature 
T* identified experimentally. In the following, we extend 
this identification to the presence of a magnetic field. The 
coupled equations to be solved are: 
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II - 0.45t||, and t± = 0.2t\^. The value of \V\ (corre- 
sponding to a given sample) is obtained by solving the 
above equations for if = and setting T*{H = 0) at 
the experimental value in zero field. These equations are 
then solved numerically for the unknowns (T* , /x) with 
given values of {n,\V\, H). In this way, to any given 
sample (corresponding to a pair of values (n, |V^|)), we 
associate a curve T*{H), from which the pseudogap clos- 
ing field _ffpg of interest is extracted by extrapolation to 
T*{Hpg) = (extrapolation is required to avoid numeri- 
cal instabilities when T*{H) approaches zero). 

A typical curve T*[H) obtained by the above proce- 
dure is shown in Fig. 1, for the values of n and \ V\ corre- 
sponding to the sample with p = 0.21 analyzed in Ref. ||. 
Note the rather weak initial dependence of T* (H) for low 
H. Typically, T*{H) decreases only by about 1% for 
values of H up to 25 T. This result is consistent , 
the experimental finding reported by several authors 
where a negligibly weak magnetic-field dependence of T* 
is noted for various compounds. 
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In these expressions, n is the site density, £,aO^) = 
—2t\\ (cos kx + cos ky) — 2t± cos k^ + 4<J| cos k^ cos ky—^ + 
HbHct (where /x is the chemical potential and a = ±1), 
2^(k) = ^+(k) -I- f-(k), 7(k) = (cosfc^; — cosfc^) for a 
d-wave gap, and each component of the wave vector k is 
integrated from — tt to +7r. To make contact with experi- 
mental data, we interpret n = 1—p where p is the doping 
valueO. Values appropriate to BSCCO are tii = 0.25 eV, 
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FIG. 1. Dependence of T*{H) on H obtained from 
Eqs. (0) and (|) with n = 0.79 and \V\ = 0.175 eV. The con- 
tinuous curve is a polynomial fit. [Temperature is expressed 
in Kelvin and magnetic field in Tesla.] 

Figure 2 shows the results of our calculation corre- 
sponding to the three samples reported in the inset of 
Fig. 4 from Ref. In that reference, the large (^ 100 
T) values of iJpg have been obtained by extrapolation of 
the data collected up to 60 T. Note that both the ex- 
perimental and the theoretical points lie on straight lines 
passing through the origin. This entails a linear pro- 
portionality between Hpg and T*{H = 0) (of the form 
2/isi7pg = akBT*{H — 0) where a is a constant), as 
pointed out in Ref. |^. [This linear proportionality has 
ben established in Ref. || not only for the three doping 
values reported in the inset of their Fig. 4, but also for 
lower doping values (albeit with larger uncertainty), as 
shown in their Fig. 4. We have also extended our cal- 
culation to the value p = 0.16, still obtaining the linear 
behavior.! 
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FIG. 2. Theoretical (squares) and experimental (circles) 
values of 2^s_ffpg/fcs vs T*{H = 0), for the three samples 
[with p = (0.20,0.21,0.22) from right to left] reported in the 
inset of Fig. 4 from Ref. 6. [Temperatures are expressed in 
Kelvin.] 

The experimental data give aoxp — 1, while our calcu- 
lation yields athco — 2.4. This difference should not be 
regarded to be too severe. First of all, recall that T* is, 
by its own nature, a non-sharply defined crossover tem- 
perature, so that there is an intrinsic uncertainty in its 
identification. Nonetheless, one may attribute the above 
discrepancy to additional physical effects (such as orbital 
frustration) not included in Eqs. (|l|) and (||). We have, 
however, verified for the continuum model (see below) 
that orbital frustration contributes only in a marginal 
way to the suppression of T* in the intermediate-coupling 
regime. In addition, one may observe that, while the ex- 
perimental value Tgxp (as obtained in Ref. ||) identifies 
directly an energy scale (which, by our interpretation, is 
associated with theJDinding energy Ef, of a pair embed- 
ded in the mediumllil), the theoretical value T^^^^^ can be 
smaller than Ei,. This statement is readily verified in the 
strong-coupling limit of the continuum model considered 
beW whereby kBT*^^^^ = {eo/2) / \a{to/ epf^ in zero 
fieldlij, ep being the Fermi energy and to the binding en- 
ergy of the associated two-body problem {E^ is expected 
to reduce to Co in this limit). From this relationship, 
one gets T*^^^ < to when Co/ep > exp(l/3) ~ 1.4. At 
weaker (intermediate) couplings, on the other hand, nu- 
merical calculation yields T^h^o — eo/2 when Eo/cf — 1- 
In this case, the experimental relation g^eHpg = kBT*^^ 
is thus equivalent to the theoretical relation gfiBHpg — 
atheo^B^thoo with athco — 2 (provided one still identifies 
El, with Eo). The difference between the experimental 
and theoretical identification of T* may thus account for 
the difference between oiexp and athco- 

In any event, the physical point to be emphasized from 
our calculation is that the energy scales associated with 
Hpg and T*{H — 0) are of the same order, as found ex- 
perimentally. The experimental result thus finds a nat- 
ural explanation within our precursor superconductivity 
scenario, whereby T*{H = 0) is associated with the bind- 



ing energy of the (quasi) preformed pairs and _ffpg with 
the Pauli pair breaking due to Zeeman splitting. 

To assess the role of orbital frustration, it is sufficient 
to analyze the simpler continuum model for which both 
Zeeman and orbital effects can be readily included. In 
this way, we avoid facing at this stage the peculiar prob- 
lems of Bloch electrons in a magnetic fieldt^rtj. 

In the continuum case, the fermionic attraction can 
be modeled by a point contact interaction, regular- 
ized (in ffhree dimensions) in terms of the scattering 
length ai?c^li3. The dimensionless interaction parame- 
ter {kpap)'^ (where the Fermi wave vector kp is related 
to the density via n = k'p/(ii:'^)) then ranges from —oo 
in weak coupling to -l-oo in strong coupling (the crossover 
region of interest being bound by [kF\ap\)^^ < 1). For 
this model, the equations identifying T* within our ap- 
proach (for an s-wave pseudogap) and corresponding to 
Eqs. (|l|) and (||) are the following: 
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where m is the fermion mass, lOc = eH/{mc) the 
cyclotron frequency, fi^ — fi — ^sHa, f{x) = 
[exp(x/fcB T*) + 1]^^ the Fermi function at temperature 
T*, and 
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N{e) being the independent-particle density of states per 
spin component and F{e) = f{e — -f /(e — /i_) — 
1. These equations contain both the Zeeman splitting 
(through the spin dependence of and the orbital ef- 
fect (through the presence of b(T*,iJ,;H) in Eq. (||) and 
of the Landau levels in Eq. (^)). These equations in- 
clude, however, the effect of the magnetic field on a dif- 
ferent footing. Equation (||) (which has been obtained 
by considering the Thouless instability in the particle- 
particle channel) bears on the eikonal approximation, 
whereby the orbital effect of the magnetic field is intro- 
duced by multiplying the real-space representation of the 
independent-particle Green's functipa in the absence of 
the field by a phase factor exp{i(/3}ll3. Equation (^), on 
the other hand, keeps the Landau orbital quantization ex- 
plicitly, otherwise the eikonal approximation would com- 
pletely wash out the magnetic orbital effect. 

It is relevant to discuss at this point the conditions for 
the validity of the eikonal approximation in the present 
context. The domain of validity of the eikonal approx- 
imation for the single-particle Green's function depends 
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on the distance between its two spatial coordinates. For 
Eq. (^) , this domain is delimited by the range I of the two- 
particle propagator in zero field, yielding ujc pc/iml) 
where pc is largest momentum scale of the independent- 
particle Green's function. In addition, by expanding up 
to second order the phase factor exp{iip} of the eikonal 
approximation in Eq. (0), one arrives at the additional 
condition ujc ^ {mP)~ . In weak coupling, pc ~ kp and 
I Ki where is the Pippard coherence length. In this 
case, one obtains lOc ^ (™Co)~"'^ f^'' t-^ic validity of the 
eikonal approximation. In strong coupling, on the other 
hand, Pq^ Ki ap Ki leading to ujc ^ Co for the validity 
of the eikonal approximation. 

The above estimates lead us to conclude that, strictly 
speaking, the eikonal approximation holds only for the 
initial part of the curve T*{H) close to zero field. Tak- 
ing tJc^Co ^ 0-1 i'^ weak and oJc/^o ^ 0.1 in strong 
coupling, one obtains correspondingly {T*{H = 0) — 
T*(H))/T*{H = 0) < 1% and < 10%. The relative 
range of validity of the eikonal approximation increases, 
therefore, when passing from weak to strong coupling. In 
practice, we shall extrapolate the results for T*{H) ob- 
tained from Eqs. (||) and outside the strict range of 
validity of the eikonal approximation to obtain the pseu- 
dogap closing field -ffpg, since we are interested in ob- 
taining a qualitative estimate of the relative importance 
of the Zeeman and orbital effects at intermediate cou- 
pling. The eikonal approximation, on the other hand, 
does not affect the equations determining T*{H) when 
the Zeeman effect only is included. 



density n/2). Note that atheo — 2.6 with the Zeeman ef- 
fect only, while atheo — 2.3 when both effects are included 
(to be compared with the value atheo — 2.4 from Fig. 2). 
This implies that the calculation with the Zeeman effect 
only is rather stable against the inclusion of the orbital 
effect. Note also from this figure that, in the small H 
region, T*{H) is somewhat sensitive to the inclusion of 
the orbital effect (albeit the dependence of T* (H) on H 
is still weak in an alisolute sense) , in agreement with the 
experimental resultB for fields up to 14 T. 

In conclusion, we have shown that consideration of the 
Zeeman splitting within the precursor superconductivity 
scenario (at intermediate coupling) reproduces the exper- 
imental relation 2/XB_ffpg — kBT*{H = 0) between the 
pseudogap closing field Hpg and the pseudogap temper- 
ature T*{H = 0) at zero fieldQ. Our results are also fully 
consistent with the occurrence oLajveak initial depen- 
dence of r*(iJ) on H at low fieldsSSI. As the sensitivity 
to a magnetic field is believed to be a crucial test for the 
physical origin of the pseudogap, our results show that 
the available experimental data are fully compatible with 
a superconducting origin of the pseudogap. 
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FIG. 3. T*{H) vs H for the continuum model with cou- 
pling value [kpap)'^ = 0.5, including the Zeeman effect only 
(squares) and the orbital effect as well (triangles). The con- 
tinuous curves are polynomial fits. 

Figure 3 shows T*{H) vs H obtained by Eqs. (||) and 
(^ for the continuum model with the value (kpap)^^ — 
0.5 in the intermediate-coupling regime, by including 
both Zeeman and orbital effects (triangles) or the Zee- 
man effect only (squares) (the temperature T* has been 
conveniently normalized in terms of the Bose-Einstein 
temperature Tbe for composite bosons of mass 2m and 
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